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Abstract 

The glycoside hydrolase family (GH) 65 is a family of inverting phosphorylases that act on ot-glucosides. A GH65 protein 
(Bsel_2816) from Bacillus selenitireducens MLS10 exhibited inorganic phosphate (Pi)-dependent hydrolysis of kojibiose at the 
rate of 0.43 s -1 . No carbohydrate acted as acceptor for the reverse phosphorolysis using B-D-glucose 1 -phosphate (BGIclP) 
as donor. During the search for a suitable acceptor, we found that Bsel_2816 possessed hydrolytic activity on BGIclP with a 
k cat of 2.8 s -1 ; moreover, such significant hydrolytic activity on sugar 1-phosphate had not been reported for any inverting 
phosphorylase. The H 2 18 0 incorporation experiment and the anomeric analysis during the hydrolysis of BGIclP revealed 
that the hydrolysis was due to the glucosyl-transferring reaction to a water molecule and not a phosphatase-type reaction. 
Glycerol was found to be the best acceptor to generate 2-O-ot-D-glucosylglycerol (GG) at the rate of 180 s _1 . Bsel_2816 
phosphorolyzed GG through sequential Bi-Bi mechanism with a k cat of 95 s -1 . We propose 2-O-a-D-glucopyranosylglycerol: 
phosphate B-D-glucosyltransferase as the systematic name and 2-O-a-D-glucosylglycerol phosphorylase as the short name 
for Bsel_2816. This is the first report describing a phosphorylase that utilizes polyols, and not carbohydrates, as suitable 
acceptor substrates. 
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Introduction 

Phosphorylases catalyze the cleavage of glycosidic bonds 
through substitution with phosphate [1,2]. These enzymes 
reversibly phosphorolyze glycosides to form corresponding mono- 
saccharide 1 -phosphates with retention or inversion of anomeric 
configuration. The sole exception is maltosyl-transferring starch 
synthase (EC 2.4.99.16) that generates maltose 1-phosphate [3]. 
All the phosphorylases reported (Table 1) exhibit strict substrate 
specificity in the phosphorolysis and regioselectivity in the reverse 
phosphorolysis. Phosphorylases have been classified in one of the 
following families in the Carbohydrate-Active Enzymes database 
(http://www.cazy.org/) based on the amino acid sequence 
similarity [4]: glycoside hydrolase families (GH) 13, 65, 94, 112, 
and 130, and glycosyltransferase families (GT) 4 and 35. Among 
these families, GH65, 94, 112, and 130 primarily comprises of 
phosphorylases that mainly catalyze the reversible phosphorolysis 
of a-glucosides, R-glucosides, B-galactosides, and fl-mannosides, 
respectively, with inversion of anomeric configuration. 

Various oligosaccharides have been synthesized by phosphor- 
ylases using the reverse phosphorolysis from the corresponding 
sugar 1-phosphate as a donor substrate and suitable carbohydrate 
acceptors [1,2,5,6]. In addition, the reversibility allows the 
practical synthesis of oligosaccharides from the readily available 



natural sugars by using a single phosphorylase [7,8], by the 
combined reaction of two phosphorylases that share the same 
sugar 1-phosphate or utilizing the compatible ones [1,9-13], or by 
the combined reaction of two phosphorylases producing different 
sugar 1 -phosphates with additional enzymes to convert the sugar 
1 -phosphates [14-16]. 

However, the utilizations of phosphorylases for the practical 
syntheses of oligosaccharides are limited by the fact that there is 
little variation among the phosphorylases, with only 22 known 
activities (Table 1). Therefore, it is of interest to identify 
phosphorylases with other substrate specificities. The following 
new phosphorylases have been reported within the past five years: 
GH13 maltosyl-transferring starch synthase (EC 2.4.99.16) [3], 
GH65 nigerose phosphorylase (EC 2.4.1.279) [17], GH65 3-0-a- 
D-glucosyl-L-rhamnose phosphorylase (EC 2.4.1.282) [18], GH94 
cellobionic acid phosphorylase [19], GH1 12 4-O-R-D-galactosyl-L- 
rhamnose phosphorylase (EC 2.4.1.247) [20], GH130 0-1,4-d- 
mannosyl-D-glucose phosphorylase (EC 2.4.1.281) [21], GH130 B- 
1 ,4-mannooligosaccharide phosphorylase (EC 2.4.1.x) [22], and 
GH130 4-0-B-D-mannosyl- l A , '-acetyl-D-glucosamine phosphorylase 
(EC 2.4.1.x) [23]. 

Bacillus selenitireducens is an anaerobic haloalkaliphilic bacterium 
that is often isolated from the bottom sediments [24]. It grows 
optimally in a salinity of 1 .5 M NaCl. The salt-dependent proteins 



PLOS ONE | www.plosone.org 



1 



January 2014 | Volume 9 | Issue 1 | e86548 



2-0-a-D-Glucosylglycerol Phosphorylase 



Table 1. 


List of known phosphorylases. 










EC 


name 


mechanism 


product 


Family 


2.4.1.1 


glycogen/starch phosphorylase 


retention 


aGlclP 


GT35 


2.4.1.7 


sucrose phosphorylase 


retention 


aGlclP 


GH13 


2.4.1.8 


maltose phosphorylase 


inversion 


PGIclP 


GH65 


2.4.1.20 


cellobiose phosphorylase 


inversion 


aGlclP 


GH94 


2.4.1.30 


1,3-P-oligoglucan phosphorylase 


inversion 


aGlclP 


Nl* 


2.4.1.31 


laminaribiose phosphorylase 


inversion 


aGlclP 


GH94 


2.4.1.49 


cellodextrin phosphorylase 


inversion 


aGlclP 


GH94 


2.4.1.64 


trehalose phosphorylase (inverting) 


inversion 


PGIclP 


GH65 


2.4.1.97 


f}-1,3-glucan phosphorylase 


inversion 


aGlclP 


Nl* 


2.4.1.211 


1 ,3-P-galactosyl-A/-acetylhexosamine phosphorylase 


inversion 


aGallP 


GH112 


2.4.1.216 


trehalose 6-phosphate phosphorylase 


inversion 


PGIclP 


GH65 


2.4.1.230 


kojibiose phosphorylase 


inversion 


PGIclP 


GH65 


2.4.1.231 


trehalose phosphorylase (retaining) 


retention 


aGlclP 


GT4 


2.4.1.247 


1 ,4-p-D-galactosyl-L-rhamnose phosphorylase 


inversion 


aGallP 


GH112 


2.4.1.279 


nigerose phosphorylase 


inversion 


PGIclP 


GH65 


2.4.1.280 


W,/v'-diacetylchitobiose phosphorylase 


inversion 


aGlcNAclP 


GH94 


2.4.1.281 


1 ,4-p-mannosyl-glucose phosphorylase 


inversion 


aManIP 


GH130 


2.4.1.282 


1 ,3-ot-D-glucosyl-L-rhamnose phosphorylase 


inversion 


PGIclP 


GH65 


2.4.1.x 


,4-mannooligosaccharide phosphorylase 


inversion 


aManIP 


GH130 


2.4.1.x 


1 ,4-p-mannosyl-/v-acetylglucosamine phosphorylase 


inversion 


aManIP 


GH130 


2.4.1.x 


cellobionic acid phosphorylase 


inversion 


aGlclP 


GH94 


2.4.99.16 


starch synthase (maltose-transferring) 


retention 


aMallP 


GH13 



*Not identified. 

doi:1 0.1 371 /journal.pone.0086548.t001 



isolated from this bacterium, such as the halophilic extracellular 
arsenate reductase that shows its optimal activity at 2.5 M NaCl, 
can be adapted to function in extreme environments [25]. The 
genomic sequence of the B. selenitireducens MLS 1 0 strain, which was 
isolated from the anoxic mud of a saline lake, has been published 
(GenBank accession No. CP001791.1). B. selenitireducens MLS 10 
possesses three genes encoding possible inverting phosphorylases 
belonging to GH65. We have revealed the activities of two of 
them: one (locus: Bsel_2056) is a maltose phosphorylase that 
produces derivatives of maltose possessing a branch at position 2 
[26], and the other (locus: Bsel_1207) is a potassium ion- 
dependent trehalose phosphorylase [27]. In this study, we 
identified the specificity of the third GH65 protein (locus 
Bsel_2816, GenBank accession No. ADI00307.1) and found that 
it is unique among the known phosphorylases. 



Materials and Methods 

Materials 

(3-D-Glucose 1 -phosphate (fiGlcl/) disodium salt was purchased 
from Tokyo Chemical Industry (Tokyo, Japan). (3-Phosphogluco- 
mutase (pTGM), PGIclP bis(cyclohexylammonium) salt, and D- 
glucose 1,6-bisphosphate (G\c\6bP) were obtained from Sigma- 
Aldrich Chemicals (St. Louis, MO). Kojibiose, glucose 6- 
phosphate dehydrogenase (G6PDH, from Leuconostoc mesenteroides), 
thio-NAD + , and Glucose C-II Test Wako were purchased from 
Wako Pure Chemicals (Osaka, Japan). All other chemicals used 
are of reagent grade. 
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Amino Acid Sequence Analysis 

Multiple alignments of the amino acid sequences of GH65 
proteins with known activities with that of Bsel_2816 were 
performed using the ClustalW version 2.1 on the DDBJ sever 
(http:// clustalw.ddbj.nig.ac.jp/index.php?lang = ja). A phyloge- 
netic tree was drawn from the results obtained with the 
neighbor-joining method using a Tree View version 1.6.6 
(http://taxonomy.zoology.gla.ac.uk/rod/rod.html). 

Cloning of bsel_2816 

A gene (GenBank ID: ADI00307.1) encoding Bsel_2816 was 
amplified by PCR from the genomic DNA of B. selenitireducens used 
as template and KOD-plus DNA polymerase (Toyobo, Osaka, 
Japan) with the following primers based on the genomic sequence: 
5'-aaaccatgggccatgaaattggagaacatc-3' as the forward primer 
containing the Ncol site (underlined) and 5 ' -ttt ctcgag gcgggacttggt- 
gatgcg-3'as the reverse primer containing the Xhol site (under- 
lined). The amplified gene was purified using the FastGene Gel/ 
PCR Extraction Kit (Nippon Genetics Co., Ltd., Tokyo, Japan), 
digested by Ncol and Xhol (New England Biolabs, Beverly, MA, 
USA), and inserted into pET28a (+) (Novagen, Madison, WI, 
USA) to add a His^-tag at the C-terminus of the recombinant 
protein. The expression plasmid was propagated in Escherichia coli 
DHSot (Toyobo), purified with the FastGene Plasmid Mini Kit 
(Nippon Genetics Co., Ltd.), and verified by sequencing (Operon 
Biotechnologies, Tokyo, Japan). 
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Site-directed Mutagenesis 

Site-directed mutagenesis was carried out by the method of 
Braman et al. [28] but using the KOD-plus DNA polymerase. The 
pET28a plasmid carrying the bsel_2816 gene was used as template 
for PCR to obtain the DNA fragment containing the mutation. 
Pairs of the mutagenic primers (5'-cagggtcctgatgcataccatgagaacg- 
3' for E475A, 5'-gttcagggtcctgatcaataccatgagaac-3' for E475Q, 
and the complementary strands) were used. The parental strands 
in the PCR products were digested with Dpnl, which specifically 
digests methylated and hemimethylated DNA, and the amplified 
DNA were retained. The nicked plasmid DNA incorporating the 
desired mutation was transformed into the E. colt DH5a. The 
expression plasmid for the mutants was purified with the FastGene 
Plasmid Mini Kit (Nippon Genetics Co., Ltd.), and verified by 
sequencing (Operon Biotechnologies). 

Expression and Purification of Recombinant Proteins 

E. coli BL21 (DE3) (Novagen) transformants harboring the 
expression plasmid for Bsel_2816 or its mutants were grown at 
37°C in 200 ml of Luria-Bertani medium (1% tryptone, 0.5% 
yeast extract, and 0.5% NaCl) containing 50 |Ig/ml of kanamycin 
up to an absorbance (Abs) of 0.6 at 600 nm. The expressions were 
induced by 0.1 mM isopropyl fS-D-thiogalactopyranoside and 
continued at 18°C for 24 h. Wet cells, harvested by centrifugation 
at 20,000 xg for 20 min, were suspended in 50 mM HEPES- 
NaOH buffer (pH 7.0) containing 500 mM NaCl (buffer A) and 
10% glycerol. The suspended cells were disrupted by sonication 
(Branson sonifier 250A; Branson Ultrasonics, Emerson Japan, 
Ltd., Kanagawa, Japan) and the supernatant, collected by 
centrifugation at 20,000 xg for 20 min, was applied to a HisTrap 
HP column (GE Healthcare, Buckinghamshire, UK) equilibrated 
with buffer A containing 1 0 mM imidazole using an AKTA prime 
system (GE Healthcare). After a wash with buffer A containing 
22 mM imidazole and the following elution using a 22—400 mM 
imidazole linear gradient in buffer A, the fractions containing the 
recombinant protein (Bsel_2816) were pooled, dialyzed against 
10 mM HEPES-NaOH buffer (pH 7.0), and concentrated 
(AMICON Ultra- 15 filter; Millipore Co., Billerica, MA, USA). 
The concentration of the wild-type, E475Q mutant, and E475A 
mutant was determined spectrophotometrically at 280 nm using a 
theoretical extinction coefficient of £= 124,110 M 'cm 1 on the 
basis of amino acid sequences [29]. The molecular masses of the 
purified proteins was estimated by SDS-PAGE (Mini-PROTEAN 
Tetra electrophoresis system; Bio-Rad Laboratories, Inc., Hercu- 
les, CA, USA) and by gel filtration (HiLoad 26/600 Superdex 
200 pg; GE Healthcare) equilibrated with 10 mM HEPES- 
NaOH buffer (pH 7.0) containing 150 mM NaCl at the flow rate 
of 0.5 ml/min and using the marker proteins for molecular weight 
determination on High Pressure Liquid Chromatography (HPLC) 
(Oriental Yeast Co., Ltd., Tokyo, Japan) as standards. 

Thin Layer Chromatography (TLC) 

TLC was performed on a TLC plate (Kieselgel 60 F254; Merck, 
Darmstadt, Germany). After spotting the samples, the plate was 
developed once with acetonitrile-water (4:1 v/v) or twice with 
acetonitrile-water (9:1 v/v). The TLC plates were then soaked in 
5% sulfuric acid-methanol solution and heated in an oven until 
the spots of the carbohydrates became sufficiently visible. 

Enzymatic Quantification of Glc and (3Glc1P 

The concentration of Glc was determined by using a glucose 
oxidase/peroxidase kit, the Glucose C-II Test Wako. Two 
volumes of the coloring reagent of the kit were added to one 



volume of the sample. The increase in the Abs 5 os was measured 
after incubation for 60 min at 37°C. The concentration of fiGlclP 
was determined by using the [SPGM/G6PDH method [30] with 
the working reagent containing 5 U/ml [SPGM, 5 U/ml G6PDH, 
4 mM thio-NAD + , 20 |Ig/ml Glcl6bP, and 10 mM MgCl 2 in 
50 mM MOPS-NaOH buffer (pH 7.0). An aliquot of sample was 
mixed with the same volume of the working reagent and the 
increase in the Abs 4 oo was measured after incubation for 60 min at 
37°C. 

Cleavage of oc-linked Glucobioses 

To examine the specificity on a-linked glucobioses, each 
reaction mixture containing each disaccharide (10 mM; kojibiose, 
nigerose, maltose, isomaltose, or trehalose), 10 mM inorganic 
phosphate (Pi), and 29 |0.M Bsel_2816 in 50 mM MES-NaOH 
buffer (pH 6.0) was incubated at 30°C for 3 h. The products were 
analyzed by TLC. The reaction rate for the cleavage of kojibiose 
was determined by quantifying the Glc liberated in the reaction 
mixtures containing various concentrations of kojibiose (1— 
10 mM) and Pi (0-10 mM) in 40 mM sodium acetate buffer 
(pH 5.5) at 30°C. An aliquot of the reaction mixture was mixed 
with the same volume of dimethyl sulfoxide to stop the enzymatic 
reaction, followed by the enzymatic quantification of Glc. The 
kinetic parameters were calculated by curve-fitting the experi- 
mental data to the theoretical equation (i) for sequential Bi Bi 
mechanism (A = kojibiose, B = Pi) using GraFit version 7.0.2 
(Erithacus Software Ltd., London, UK). 

v = k CAt [E] 0 [A][B}/(K iA K mB +K mA [B,}+K mB [A] + [A][B]) (i) 

Acceptor specificity analysis. To determine the acceptor 
specificity of Bsel_2816 in the reverse phosphorolysis, each 
reaction mixture containing 10 mM each acceptor and 10 mM 
PGlclPin 50 mM MES-NaOH buffer (pH 6.0) was reacted with 
Bsel_2816 at 30°C. After the reaction, the products were analyzed 
by TLC. 

Preparation and Structural Determination of the Product 
from Glycerol 

A reaction mixture (5 ml) containing 500 mM (3GlclP (disodi- 
um salt buffered at pH 6.0 with HC1), 1 M glycerol, and 0.25 mg/ 
ml Bsel_2816 was incubated at 30°C for 18 h. After terminating 
the reaction by boiling, an aliquot (3 ml) was deionized with 
Amberlite MB-3 (Organo, Tokyo, Japan), applied on an Amberlite 
IRA400 column (OH" type, 40 ml), and eluted with H 2 0 to 
remove the reducing sugars. The fractions containing the product 
were collected, concentrated to 9 ml, and loaded on a Toyopearl 
HW40F column (50 mm internal diameter x 800 mm; Tosoh, 
Tokyo, Japan) equilibrated with distilled water. The product was 
eluted by distilled water at a flow rate of 5 ml/ min. The fractions 
containing the product were collected, concentrated, and lyoph- 
ilized to be syrup (298 mg). The one-dimensional ('H and L5 C) 
and two-dimensional [double-quantum-filtered correlation spec- 
troscopy (DQF-COSY), heteronuclear single-quantum coherence 
(HSQC), and heteronuclear multiple-bond correlation (HMBC)] 
NMR spectra of the product were recorded in D 2 0 with 2-methyl- 
2-propanol as an internal standard at 298 K using the Bruker 
Avance 500 spectrometer (Bruker Biospin, Rheinstetten, Ger- 
many). The product was determined to be 2- O-ot-D-glucosylgly- 
cerol (GG). 
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Measurement of Enzyme Activity 

The phosphorolytic activity of GG was determined by 
quantifying (5Glc 1 P produced during the enzymatic reaction from 
10 mM GG and 10 mM Pi in 40 mM HEPES-NaOH buffer 
(pH 8.0) at 30°C An aliquot of the reaction mixture was mixed 
with the same volume of dimethyl sulfoxide to stop the enzymatic 
reaction, followed by the enzymatic quantification of PGlclP. The 
hydrolytic activity on PGlc 1 P was determined by measuring the 
increase in Glc released during the enzymatic reaction from 
10 mM PGlclP in 40 mM MES-NaOH buffer (pH 6.0) or 
40 mM HEPES-NaOH buffer (pH 7.5) at 30°C. An aliquot of the 
reaction mixture was mixed with the same volume of dimethyl 
sulfoxide to stop the enzymatic reaction, followed by the 
enzymatic quantification of Glc. The synthetic activity of GG 
was determined by measuring the increase in Pi in the reaction 
mixture containing 10 mM PGlclP and 10 mM glycerol in 
40 mM HEPES-NaOH buffer (pH 7.5) at 30°C by using the 
method of Lowry and Lopez [31]. 

Temperature and pH Profile 

The effects of pH on the phosphorolytic, hydrolytic, and 
synthetic activities of 212 nM Bsel_2816 were measured under the 
standard conditions by substituting each 40 mM reaction buffer 
into the following 40 mM buffers: sodium citrate (pH 3.0—5.5), 
Bis(2-hydroxyethyl)iminotris(hydroxymethyl)methane-HCl 
(pH 5.5-7.0), HEPES-NaOH (pH 7.0-8.5), and glycine-NaOH 
(pH 8.5-10.5). The thermal and pH stabilities were evaluated by 
measuring the residual synthetic activity under the standard 
conditions after the incubation of Bsel_2816 (7 |0,M and 11 |0M, 
respectively) at temperatures 30-90°C for 15 min in 10 mM 
HEPES-NaOH buffer (pH 7.0) containing 0.1% bovine serum 
albumin and at the various pH values at 4°C for 24 h, respectively. 

The Kinetic Parameters of Bsel_2816 

A kinetic analysis of the phosphorolytic reaction of GG was 
performed using the continuous PGlclP assay. The reaction 
mixtures (150 (0.1) were prepared in the wells of a 96-well microtiter 
plate with various concentrations of GG (1.0-20 mM) and Pi (0.5- 
10 mM), and 5.2 nM Bsel_2816 in 40 mM HEPES-NaOH buffer 
(pH 8.0) containing 2.5 U/ml pTGM, 2.5 U/ml G6PDH, 2 mM 
thio-NAD + , 5 |ig/ml Glcl6bP, and 5 mM MgCl 2 . The reaction 
was carried out in a temperature-controlled microplate reader 
(Multiskan GO, ThermoFisher Scientific, Waltham, MA) at 30°C, 
and the Abs4oo was continuously monitored at 30-s intervals 
without stopping the reaction. The kinetic parameters were 
calculated by curve-fitting the experimental data to the theoretical 
equation (i) for sequential Bi Bi mechanism (A = GG, B = Pi) using 
the GraFit version 7.0.2. 

The initial velocities of the hydrolytic reaction on PGlclPwere 
determined under the standard conditions with 5.2 nM Bsel_2816 
and various concentrations of PGlclP (0.1-10 mM). The kinetic 
analysis of the synthetic reaction with glycerol as acceptor was 
carried out under the standard conditions (5.2 nM Bsel_2816) and 
various concentrations of glycerol (1.0-100 mM) or PGlclP 
(ranging from 0.1 to 10 mM) with 10 mM or 100 mM of each 
opposite substrate, respectively. The kinetic parameters were 
calculated by curve-fitting the experimental data with the 
Michaelis-Menten equation (ii) using GraFit version 7.0.2. 

v = A: cat [E] 0 [S]/( J S: m + [S]) (ii) 



The inhibition of the hydrolysis of PGlclP by glycerol was 
examined in the reaction mixture containing 10 mM PGlclPand 
various concentrations of glycerol (0-100 mM) in 40 mM 
HEPES-NaOH buffer (pH 7.5) and the increase in the concen- 
tration of Glc was quantified. The rate of the release of Glc was 
regressed with the concentration of glycerol to determine the 
dissociation constant of glycerol using the equation (iii) with GraFit 
version 7.0.2. 

v = v 0 -K d /(K d + [glycerol]) (iii) 



Position of Hydrolysis in PGlclP 

PGlclP bis(cyclohexylammonium) salt (46 mg, 0.10 mmol) was 
dissolved in 1 ml of H 2 0 and applied to a Dowex 50 column (H + 
form, Wako Pure Chemicals) to remove the counter- ion (cyclo- 
hexylamine). An appropriate amount of 1 M NaOH was added to 
the decationized solution containing free PGlclP (5 ml) adjusting 
pH to 6.0, followed by concentration to 1 ml by a centrifuge 
evaporator to obtain 100 mM PGlclP-Na buffer (pH 6.0). An 
aliquot of the 100 mM PGlclP-Na buffer (20 |0,1) was mixed with 
70 |ll of H 2 18 0 (97%, Sigma-Aldrich) or unlabeled H 2 0 as a 
control. The hydrolysis of PGlc 1 P was started by adding 1 0 |ll of 
29 uM Bsel_2816 in 4 mM MES-NaOH buffer (pH 6.0) to the 
solution. After the incubation at 30°C for 4 h, PGlclP was 
completely hydrolyzed into Glc and Pi. The solution was 
lyophilized, dissolved in 50% methanol, and subjected to negative 
and positive ESI-MS analysis by using a Bruker APEX II 70e 
Fourier transformed cyclotron resonance mass spectrometer 
(Bruker Daltonics, Billerica, MA, USA). 

Specificity on the Anomeric Configuration 

The anomeric configuration of the glucose generated during the 
hydrolysis of PGlclP was analyzed by HPLC [32], using an HPLC 
system (Prominence, Shimadzu, Kyoto, Japan) equipped with a 
refractive index detector, RID 1 OA (Shimadzu). The hydrolysis of 
PGlclPwas performed in the reaction mixture containing 50 mM 
PGlclP and 13.8 (iM Bsel_2816 prepared in 20 mM MES- 
NaOH buffer (pH 6.0) at 30°C. After the reaction for 1.5 and 
25 min, aliquots (10 J0.1 each) of the reaction mixtures were 
immediately injected into the HPLC and separated with a TSK 
amide-80 column (4.6 mm internal diameter x 250 mm; Tosoh, 
Tokyo, Japan) by isocratic mode with 80% acetonitrile at a flow 
rate of 1.5 ml/min. 

Results 

Gene Cloning, Expression, and Purification of 
Recombinant Bsel_2816 

The phylogenetic tree analysis showed that Bsel_2816 is 
distandy related with the kojibiose phosphorylase and nigerose 
phosphorylase (Fig. 1). In this study, the bsel_2816 gene cloned 
from the genomic DNA of B. selenitireducens MLS 1 0 was expressed 
as a Hise-tag fusion protein in E. coli BL2 1 (DE3) to investigate its 
enzymatic properties. Purified protein (approximately 14 mg) was 
obtained from the cell lysate of 200 ml culture. The purified 
Bsel_2816 migrated in SDS-PAGE as a single protein band with 
an estimated molecular mass of approximately 87 kDa, in 
agreement with the theoretical molecular weight of 86,259. 
However, this native molecular mass was determined to be 
174 kDa by gel filtration, suggesting that Bsel_2816 was a 
homodimeric enzyme. 
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Figure 1. Phylogenetic tree of the GH65 enzymes characterized. Sequences used are: 1, Bsel_2816 of B. selenitireducens MLS10 (GenBank ID: 
ADI00307.1); 2, AII1058 of Nostoc sp. PCC 7120; 3, AII4989 of Wosfoc sp. PCC 7120; 4, Csac_0439 of Caldicellulosiruptor saccharolyticus DSM 8903 
(ABP66077.1); 5, KojP of Thermoanaerobacter brockii ATCC 35047 (AAE30762.1); 6, Bsel_2056 of B. selenitireducens MLS10 (ADH99560.1); 7, MapA of 
Paenibacillus sp. SH-55 (BAD97810.1); 8, maltose phosphorylase of Lactobacillus sanfranciscensis DSM 20451 (CAA1 1905.1); 9, MPase of Bacillus sp. RK- 
1 (BAC54904.1); 10, EF0957 of Enterococcus faecalis V583 (AAO80764.1); 11, LBA1870 of Lactobacillus acidophilus NCFM (AAV43670.1); 12, maltose 
phosphorylase from Lactobacillus brevis ATCC 8287 (Q7SIE1); 13, TreP of T. brockii ATCC 35047 (AAE18727.1); 14, Bsel_1207 of B. selenitireducens 
MLS10 (ADH98720.1); 15, TPase of Geobacillus stearothermophilus SK-1 (BAC20640.1); 16, Cphy_1874 of Clostridium phytofermentans ISDg 
(ABX42243.1); 17, Cphy_1 01 9 of C. phytofermentans ISDg (ABX41 399.1); 18, TrePP of Lactococcus lactis subsp. lactis 111403 (AAK04526.1); 19, Atml of 
Metarhizium acridum CQMA102 (ABB51 158.1); 20, TreA of Aspergillus nidulans FGSCA4 (EAA66407.1 ); 21, Ath1 of Saccharomyces cerevisiae S288c 
(CAA58961.1); 22, Atc1 of Candida albicans (AAV05390.1). Proteins from 6. selenitireducens MLS10 are underlined. 
doi:10.1371/journal.pone.0086548.g001 



Examination of the Phosphorolytic Activity on the 
Substrates of the known GH65 Phosphorylases 

The phosphorolytic activity of Bsel_2816 on ot-linked gluco- 
bioses (i.e., kojibiose, nigerose, maltose, isomaltose, and trehalose) 
including the known substrates of the GH65 enzymes was 
examined. Bsel_2816 cleaved kojibiose at a slow rate only in the 
presence of Pi. To our surprise, Glc was detected as the sole 
product and (3GlclPwas not detected at all. The optimum pH for 
the cleavage of kojibiose was determined to be 5.5. The 
measurement of the reaction rate with various concentrations of 
kojibiose and Pi revealed that the reaction proceeded with a 
sequential Bi Bi mechanism, similarly to other inverting phos- 
phorylases, with the following kinetic parameters: K^j^, 8.9 mM; 
K A , 41 mM; iT mB , 0.19 mM; k vM 0.43 s" 1 , where A and B 
represent kojibiose and Pi, respectively (Fig. 2). The reaction 
mechanism suggests that Pi was involved in the cleavage of 
kojibiose even though no fJGlclPwas detected in the products. 



Hydrolytic Activity on p"Glc1P 

To determine the substrates of Bsel_2816 as a GH65 
phosphorylase, the reverse phosphorolysis was examined with 
various sugars [hexoses (D-glucose, D-mannose, D-allose, D- 
galactose, D-altrose, D-talose, and D-fructose), pentoses (D-xylose, 
D-lyxose, D-ribose, D-arabinose, and L-arabinose), aminosugars (d- 
glucosamine, D-galactosamine, jV-acetyl-D-glucosamine, jV-acetyl- 
D-mannosamine, jV-acetyl-D-galactosamine), uronic acids (D-glu- 
curonic acid, and D-galacturonic acid), deoxyhexoses (1,5-anhy- 
dro-D-glucitol, 2-deoxy-D-glucose, and L-rhamnose), derivatives of 
D-glucose (oc-D-glucose 1-phosphate, D-glucose 6-phosphate, 3-0- 
methyl-D-glucose, methyl-oc-D-glucoside, and methyl-fi-D-gluco- 
side), and disaccharides (trehalose, kojibiose, nigerose, maltose, 
isomaltose, sucrose, sophorose, laminaribiose, cellobiose, gentio- 
biose, lactose, lactulose, melibiose, xylobiose, P-l,4-mannobiose, 
and jV,jV'-diacetylchitobiose)] as acceptors and fJGlclP as the 
donor. No transfer product was detected on TLC with any sugar 
acceptors that were examined. However, the consumption of 
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Figure 2. Double-reciprocal plot of the Pi-dependent hydroly- 
sis of kojibiose by Bsel 2816 at various concentrations of Pi. 

Open circles, 1 mM Pi; closed circles, 2 mM Pi; open squares, 3 mM Pi; 
closed squares, 5 mM Pi; and open triangles, 10 mM Pi. The kinetic 
parameters were k cat = 0.43 + 0.02 (s -1 ), K mA = 8.9±0.6 (mM), 
K mB = 0.19±0.12 (mM), and K iA = 41±29 (mM). The values were 
determined by regressing the data with the following equation using 
the Grafit Version 7.02: v = WE]o[A][B]/(K iA K mB +/WB]+K mB [A]+[A][B]). 
doi:1 0.1 371 /journal.pone.0086548.g002 

fSGlclPand the concomitant generation of Glc were detected with 
all the acceptors. The concomitant generation of Glc and Pi was 
detected during the reaction in the absence of the acceptor 
molecule, clearly indicating that Bsel_2816 hydrolyzed PGlclP. 
The rates of the hydrolysis of 1 0 mM |3Glc 1 P in the presence and 
in the absence of 1 0 mM Glc were identical. The optimum pH for 
the hydrolysis was determined to be 6.0. The kinetic parameters of 
the hydrolysis of PGlclP at pH 6.0 and pH 7.5 were determined 
to be as follows: K m , 0.041 and 0.16 mM; k cM , 2.8 and 1.7 s~\ 
respectively. 

Position of the Hydrolysis in pGlclP 

The enzymatic hydrolysis of phosphate esters is generally 
catalyzed by phosphatases that cleave the PO bond of the 
phosphate ester [33]. To investigate the mechanism of the 
hydrolysis of PGlclP by Bsel_2816, the incorporation of an 
oxygen atom was examined in O-labled water. In negative MS 
(Fig. SI A) of the control reaction, the ion signals at 194.95 ml 'z 
corresponding to [2H 3 P0 4 - H]~ (calcd 194.947) and 277.03 ml 
z corresponding to [H 3 P0 4 + Glc — H] (calcd 277.033) were 
detected. In the H 2 ls O reaction, no signal of [2H 3 P0 4 — H 
( ls O)]" (calcd 196.951) and [2H 3 P0 4 - H (2 1!i O)]" (calcd 
198.955) was detected, clearly indicatrng that 18 0 was not 
incorporated in the Pi during the hydrolysis. On the other hand, 
the signal at 277.03 ml z decreased and a new signal at 279.04 ml z 
corresponding to [H 3 P04 + Glc — H ( 18 0)]~ appeared, indicating 
that O was incorporated into either Pi or Glc. In positive MS 
(Fig. SIB), a signal at 203.05 ml 'z corresponding to [Glc+Na] + 
(calcd 203.053) accompanied by a signal at 205.06 ml z corre- 
sponding to [Glc+Na( ls O)] + indicated that ls O was incorporated 
into Glc. It is concluded that the anomeric CO bond was cleaved 
during the hydrolysis of PGlclP, indicating that the hydrolysis was 
a glycosidase-type reaction and not a phosphatase-type reaction 
(Fig. 3). 




Figure 3. Linkages cleaved in the hydrolysis of pGlclP. A: 

Cleavage positions by the action of phosphatases and glycosidases. B: 
Results of the enzymatic hydrolysis of fjGlclP by Bsel_2816 in H 2 ,s O. 
doi:10.1371/journal.pone.0086548.g003 

The Anomeric Configuration of the Hydrolytic Product of 
PGIdP 

The anomeric configuration of Glc generated through the 
hydrolysis of PGlclP by Bsel_2816 was determined by HPLC 
(Fig. 4). The oc-anomer of Glc was predominant at the initiation of 
PGlclP hydrolysis, indicating that Bsel_2816 catalyzes the 
hydrolysis of PGlclP with the inversion of the anomeric 
configuration. This result suggests that the hydrolysis was not 
due to a phosphatase-type reaction, in which the anomeric 
configuration of Glc should be retained during the hydrolysis. In 
addition, the anomeric inversion suggests that water molecules 
play a role as acceptors in the reverse phosphorolysis causing the 
hydrolysis. 

Determination of the Specificity of Bsel_2816 

Because no acceptor was found within sugars, various alcohols 
(1-butanol, 2-butanol, cyclohexanol, benzylalcohol, hydroquinone, 
and phenol) and polyols (ethylene glycol, glycerol, erythritol, 
xylitol D-arabinitol, L-arabinitol, and D-glucitol) were used to 
identify the suitable acceptor. Only the hydrolysis of PGlclP was 
detected with all the alcohols examined (data not shown). Among 




0 5 10 15 

Retention time (min) 



Figure 4. The anomeric specificity during the hydrolysis of 
PGlclP. 

doi:10.1371/journal.pone.0086548.g004 
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polyols, the consumption of PGlcl-P was obviously faster with 
glycerol than with others (Fig. 5). A spot exhibiting slightly 
different Rvalue with Glc was detected on TLC. The synthetic 
reaction gave a single product from glycerol. The NMR spectra 
are shown in Fig. S2 and the assignment of the signals are 
summarized in Table 2. The Ji ;2 value of the Glc residue (3.9 Hz) 
indicates an a-D-glucoside. The HMBC cross signals indicate the 
glycosyl linkage with the position 2 of glycerol. Thus, the product 
was identified to be GG. 

In the presence of glycerol, the synthesis of GG was much faster 
than the hydrolysis of PGlclP (Fig. 6). In the presence of 10 mM 
glycerol, the hydrolysis rate was only 0.63% of that of the synthesis 
of GG. The hydrolysis was inhibited by glycerol (Fig. 6) with the 
value of 7.6 mM. In addition, Bsel_2816 exhibited a weak 
reverse phosphorolysis activity on ethylene glycol, yielding two 
products (Glc and unknown) on TLC (Fig. 5). We confirmed the 
phosphorolysis of GG by the enzyme. During the phosphorolysis, 
the major product was (3GlclP, but a trace amount of Glc was 
detectable (Fig. 7). Herein, we propose 2-0-<x-D-glucopyronosyl- 
glycerol: phosphate P-D-glucosyltransferase as the systematic name 
and 2-O-a-D-glucopyronosylglycerol phosphorylase (GGP) as the 
short name for Bsel_2816. 

Basic Properties of Bsel_2816 

Bsel_2816 was stable up to 35°C during 15-min incubation, and 
it was stable at pH between 5.5 — 10.5 at 4°C for 24 h (data not 
shown). The optimum pHs for the phosphorolytic, synthetic, and 
hydrolytic reactions were pH 8.0, pH 7.5, and pH 6.0, respec- 
tively (data not shown). 

Kinetic Analyses of Bsel_2816 

The phosphorolytic reaction of GG was investigated at pH 8.0 
at 30°C. Double reciprocal plots of the initial velocities against 
various initial concentrations of GG and P; presented a series of 
lines intersecting at a point (Fig. 8), indicating that the 
phosphorolysis follows a sequential Bi Bi mechanism, as reported 
with other inverting phosphorylases [1,17,19,20,26,34—39]. The 
kinetic parameters were calculated to be K^a, 1.1 mM; K^, 
18 mM; K mB , 0.57 mM; and k cM 95 s~\ where A and B are GG 
and Pi, respectively. The parameters for GG are in the same 
ranges of those of the other inverting phosphorylases 





S1 S21 2 3 4 

Figure 5. Screening for the suitable acceptor within polyols in 
the reverse phosphorolysis by Bsel 281 6. The acceptors used are: 
1, ethylene glycol; 2, glycerol; 3, erythritol; 4, no acceptor. SI and S2, Glc 
and as fiGlcl P as standards, respectively. The TLC plate was developed 
twice with acetonitrile-water (9:1 v/v). 
doi:1 0.1 371 /journal.pone.0086548.g005 



Table 2. Assignments of 
glucosylglycerol. 


H and 13 C NMR 


signals of 2-O-oc-d- 




Residue Position 


13 C-NMR 


'H-NMR 








5 (ppm) 


5 (ppm) 


pattern 


-MHzl 


glucose 1 


99.55 


5.11 


d 


3.9 


2 


73.21 


3.54 


dd 


3.9, 9.9 


3 


74.60 


3.73 


m 




4 


71.26 


3.40 


t 


9.5 


5 


73.67 


3.83 


m 




6 


62.18 


3.83 


m 




6' 




3.75 


m 




glycerol 1 


63.10 


3.74 


m 




r 




3.74 


m 




2 


80.49 


3.80 


m 




3 


62.04 


3.76 


dd 


4.3, 12.1 


3' 




3.70 


dd 


5.2, 12.1 



Bold: HMBC cross signal is observed with the anomeric position of glucose. 
doi:1 0.1 371 /joumal.pone.0086548.t002 



200 r- 1 1 1 1 1 1 1 1 1 J 2 




[glycerol] (mM) 

Figure 6. The reaction rates of the generations of GG and Glc 
from 10 mM fiGlcl ^ and various concentrations of glycerol. 

Open circles, d[GG]/dt; closed circles, d[Glc]/dt. 
doi:10.1371/journal.pone.0086548.g006 

[1,17,19,20,26,34-39], indicating that GG is the true substrate 
of Bsel_2816. 

The kinetic parameters of the substrates in the synthetic 
reaction (glycerol and PGlcl-f) were determined at pH 7.5. The 
K m and k cM values with glycerol were 5.7 mM and 163 s , 
respectively. The ii m value of glycerol agreed with the K A value 
determined by the inhibition of hydrolysis of PGlclP (7.6 mM). 
The K m and k cat values of PGlclP were 0.45 mM and 180 s , 
respectively. 

Phosphorolytic and Hydrolytic Activities of Mutant 
Bsel_2816 

Alignment of the amino acid sequences of GH65 phosphory- 
lases suggests that E475 of Bsel_2816 is the general acid (for 
phosphorolysis)/base (for reverse phosphorolysis) catalyst [40,41]. 
The phosphorolytic and hydrolytic activities of two mutants of 
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Time (min) 



Figure 7. Time course of the increase of pGlcl P and Glc during 
the phosphorolysis of GG. Open circles, [BGIdP]; closed circles, 
[Glc]. 

doi:1 0.1 371 /journal.pone.0086548.g007 

Bsel_2816 at E475, E475A and E475Q, were less than 0.01% 
compared with those of the wild-type Bsel_2816, suggesting that 
both reactions were catalyzed by the same catalytic residue E475. 

Discussion 

We have determined that Bsel_2816 is a GGP. The reactions 
catalyzed by GG are summarized in Fig. 9. GGP is a unique 
phosphorylase in the following two characteristics. 

First, any other phosphorylase reported exhibits its highest 
phosphorolytic activity on the glycosyl linkages bound with sugar 
residues (Table 1). Contrarily, GGP phosphorolyzes the glucosyl 
linkage bound with glycerol. Even though GG can be produced by 
other phosphorylases such as sucrose phosphorylase [42,43], it is 
one of their minor substrates. The phosphorylases generally 
exhibit K m values 1-10 mil with their main acceptors. Structural 
analyses of various phosphorylases have revealed that these 
enzymes recognize their acceptor molecules through their 
pyranose rings (or furanose ring in the case of sucrose phosphor- 
ylase) by a number of hydrogen bond interactions to obtain their 
strict acceptor specificities with the K m values. The K m value of 
GGP with glycerol was 5.7 mM, similarly to that of other 
phosphorylases, even though glycerol is a rather small molecule 
and does not adopt any rigid conformation. 

Second, the other feature of GGP is the hydrolysis of PGlc 1 P. 
Other inverting phosphorylases do not hydrolyze sugar 1- 
phosphates [2] and this is the first report describing the significant 
hydrolysis by an inverting phosphorylase. 

To investigate the mechanism of the hydrolysis, we determined 
the linkage cleaved in the hydrolysis, the anomeric configuration 
of the hydrolytic product, and the catalytic residue of the enzyme. 
Our results prove that the hydrolysis of PGlc IP occurs through the 
same catalytic mechanism as the reverse phosphorolysis, not 
through a phosphatase-type reaction. The hydrolysis is explainable 
as a glucosyl transfer reaction from PGlclP to a water molecule 
with the anomeric inversion. The phosphate-dependent hydrolysis 
of kojibiose is not considered to be a direct hydrolysis. It is 
explained as follows: GGP phosphorolyzes kojibiose at a very slow 
rate to form PGlc 1 P that is immediately hydrolyzed to generate 
Glc and Pi by the same enzyme. In this mechanism, the ot-anomer 
of Glc is considered to be generated during the hydrolysis, 
indicating a retaining hydrolysis. Unfortunately, the anomeric 
analysis during the hydrolysis of kojibiose was unsuccessful because 




-0.4 -0.2 0 0.2 0.4 0.6 0.8 1 



1/[GG] (mM" 1 ) 

Figure 8. Double-reciprocal plot of the phosphorolysis of GG 
by Bsel 2816 at various concentrations of Pi. Open circles, 
0.5 mM Pi; closed circles, 1 mM Pi; open squares, 2 mM Pi; closed 
squares, 3 mM Pi; open triangles, 5 mM Pi; and closed triangles, 10 mM 
Pi. The kinetic parameters were /c cat = 95±2 (s _1 ), K mA =1.1±0.1 (mM), 
K mB = 0.57±0.07 (mM), and K iA =18±3 (mM). The values were deter- 
mined by regressing the data with the following equation using Grafit 
Version 7.02: v = /t cat [E] 0 [A][B]/(K iA /< mB +/C mA [B]+/( m B[A]+[A][B]). 
doi:10.1371/journal.pone.0086548.g008 

of its very slow reaction. On the other hand, the hydrolytic 
reaction was much slower than the phosphorolysis and the 
synthesis of GG, suggesting that the hydrolytic activity may be 
negligible in its in vivo role. 

Recently, indication of faint hydrolytic activity on a-mannnose 
1 -phosphate [oMa.n\ P) by UhgbMP belonging to GH130, a 0-1,4- 
mannooligosaccharide phosphorylase from unknown human gut 
bacterium, has been reported [44] . In the case, the enzyme at high 
concentration produced (5- 1 ,4-mannooligosaccharides with high 
degree of polymerization only from ocManlP, suggesting that the 
enzyme initially hydrolyzed ocManlP to generate mannose. The 
resultant mannose was not detected because it readily acted as the 
acceptor molecule to form P-l,4-mannobiose, which was a better 
acceptor than mannose to generate the polymer. This enzyme may 
share the same hydrolytic mechanism with GGP. 

GG has been reported as a compatible solute for the control of 
the cytoplasmic osmotic pressure of cyanobacteria [45,46]. 
Because B. selenitireducens is a halophilic bacterium, which is often 
isolated from soil with high concentration of salts [24] , GG may 
possess a similar role in this bacterium. The biosynthetic pathway 
of the GG of cyanobacteria is explained with the concomitant 
reactions of glucosylglycerol-phosphate synthase (EC 2.4.1.213) 
[47] and glucosylglycerol phosphatase (EC 3.1.3.69) [48]. Other 
possible enzymes that form oc-glucosyl linkage at the position 2 of 
the derivatives of glycerol are glucosyl-3-phosphoglycerate syn- 
thase (EC 2.4.1.266) [49] and glucosylglycerate synthase (EC 
2.4.1.268) [50]. However, no homologous genes encoding these 
enzymes were found in the genomic sequence of B. selenitireducens 
MLS 10 (GenBank ID: CP000896). The role of GGP is unclear but 
it may be involved in the biosynthesis of GG as well as in its 
cleavage. B. selenitireducens possesses potassium ion-dependent 
trehalose phosphorylase, even though the strain does not have 
either intracellular or extracellular biosynthetic pathways of 
trehalose [27]. B. selenitireducens has one more GH65 enzyme, 
maltose phosphorylase [26], which is a possible candidate to 
produce PGlc IP. These three enzymes may relate with the 
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Figure 9. The reactions catalyzed by GGP. A, The reaction mechanism of GGP; B, Detectable reactions catalyzed by GGP. 
doi:10.1371/journal.pone.0086548.g009 



production of GG and trehalose in the cytoplasm of B. 
selenitireducens as solutes for the control of the cytoplasmic osmotic 
pressure. 

GG is often isolated from the traditional Japanese food, 
fermented with "koji" (Aspergillus otyzae) [51]. It is valued for its 
application as a moisturizing agent and as external reagent 
protecting proteins against denaturation [43,52]. The compound 
is commercially produced from maltose and glycerol using the 
transglycosylation activity of ot-glucosidase as a mixture of ot- 
glucosylglycerols [53], or from sucrose and glycerol by the action 
of sucrose phosphorylase [43,54]. GGP is expected to be an 
alternative catalyst for the production of GG. 

Supporting Information 

Figure SI ESI-MS spectra of the hydrolytic products 

from pGlclP. A, negative mode; B, positive mode. 

(PDF) 

Figure S2 NMR spectra of the product from glycerol. 

The product was identified to be 2-0-a-D-glucopyranosylglycerol. 



The spectra were taken in D z O at 298 K with 2-methyl-2- 
propanol as an internal standard using Bruker Avance 500 
spectrometer. (A) 'rl-NMR spectrum; (B) 13 C-NMR (proton 
decupling and DEPT135); (C) DQF-COSY; (D) HSQC; (E) 
HMBC. 
(PDF) 
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